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Abstract 

We develop an analytic approach which allows us to study the behaviour 
of spin models with competing interactions and p-fold spin anisotropy, D, in 
the limit where the pinning potential which results from D is large. This 
is an expansion in inverse spin anisotropy which must be carried out to all 
orders where necessary. Interesting behaviour occurs near where the boundary 
between different ground states is infinitely degenerate for infinite D. Here as 
D decreases and the spins are allowed to soften, we are able to demonstrate 
the existence of several different behaviours ranging from a single first-order 
boundary to infinite series of commensurate phases. The method is illustrated 
by considering the soft chiral clock model and the soft clock model with first- 
and second- neighbour competing interactions. In the latter case the results 
are strongly dependent on the value of p. 
PACS 05.50.+q; 64.60.Cn; 75.10.Hk 
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I. INTRODUCTION 



There are many examples of long-period phases in nature. These include the ferrimag- 
netic phases of the rare earths, long-period atomic ordering in binary alloys and polytypism, 
the possibility of many different forms of long wavelength structural order in some minerals. 

The underlying physical mechanism responsible for the formation of long-period struc- 
tures is often the existence of competing interactions. Therefore it is of interest to understand 
the properties of model systems with such competition. Perhaps the simplest of these is the 
ANNNI model, a ferromagnetic Ising model with second-neighbour antiferromagnetic in- 
teractions along one lattice axis. Despite its simplicity the ANNNI model has a very rich 
phase structure with infinite sequences of commensurate and incommensurate phases at 
finite temperaturesSi. 

If continuous spins are considered similar structures can occur even in the ground state. 
Banerjea and Taylor! performed numerical work on the chiral X-Y model with p = 2-fold 
spin anisotropy D and showed that long-period phases are stable. Chou and Griffiths! 
later proved that an infinite number of commensurate phases appear as ground states for 
p > 3. Numerical work has also shown that the X-Y model with first- and second-neighbour 
competing interactions and p = 6-fold spin anisotropy has a highly complicated ground-state 
phase diagram!. 

Note that a common feature of these models is the spin anisotropy D. As D is increased 
from zero to infinity the spins are confined to increasingly deep potential wells and the model 
crosses over from a continuous to a discrete spin limit. For example the X-Y model with 
p-fold spin anisotropy becomes a p-state clock model for D = oo. 

For infinite spin anisotropy the ground state typically comprises a few short-period 
phases. The boundaries between the phases can either correspond to a first-order tran- 
sition where only the neighbouring phases are stable or a multiphase point at which an 
infinite number of phases are degenerate!. As D decreases from infinity the states at the 
multiphase point can either remain degenerate with the point becoming a multiphase line 
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or their energy can be differentiated as the spins soften. In the latter case some or all of 
the degenerate phases may become ground states in their own right and typically a fan of 
phases springs from the multiphase point. Indeed we shall see that in many respects D^ 1 
behaves in a way akin to temperature with the spin softening playing the part of entropic 
fluctuations. 

The large D region of the phase diagram is difficult to explore numerically because 
the phase sequences can be very complicated and the widths of the stable phases small. 
Therefore our aim in this paper is to describe an analytic technique which is useful in 
helping to understand the ground state of models with competing interactions for large spin 
anisotropy D. This is an expansion in 1/D taken to all orders where necessary. A short 
paper has summarised some of the results of the calculations^. Here our aim is to describe 
the technical details of the expansion. 

In Section |I| of the paper we explain the approach in some detail for the chiral X-Y model 
with p-fold spin anisotropy. In [11 A| the model is defined and our notation introduced. The 
energy differences which are central to the argument are defined in |11 ij| and their dependence 
on the deviation of the spins from their positions at D = oo is calculated. [11 C| descibes the 
convenient labelling of the spin states which allows a calculation of the energy differences to 
leading order presented in [II D| . From these we are able to show, in agreement with Chou 
and Griffiths! that all possible phases are stable near the multiphase point for p > 3. In 



Section [II E| the widths of the long period phases are calculated. 

The calculation is repeated in Section |I| for the X-Y model with first- and second- 
neighbour competing interactions and p-fold spin anisotropy. This is an involved calculation 
because of the existence of the second-neighbour interactions. Previous numerical results for 
this model for p = 6 were unable to proble the large D limiti. We find that the behaviour 
near the multiphase points is complicated and highly dependent on p. 

The results are summarised and discussed in Section [TV]. 
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II. THE CHIRAL XY MODEL 



A. Definitions and notation 

The approach is most easily explained by considering its application to the chiral XY 
model with p-fold spin anisotropy. This is described by the Hamiltonian 

ft = £{- J c °s(^-i - Oi + A) - D(cos p9i - l)/p 2 } (1) 

% 

where the 9i are angular variables which can take values between and 2n lying on the sites 
i of a one-dimensional lattice. 

Note that the Hamiltonian (fl]) is invariant under the transformation 

A -> A' = A + 2vrm/p, (2) 

for any integer m given the reidentification 

df, => Q\ = (9i + 2mni/p) (3) 

and therefore we may restrict our attention to < A < 2Tr/p. Moreover the system is 
invariant under 

A A' = 2ir/p - A, 

e j ^e' j = (-e j + 2nj/p). (4) 

Thus the phase boundaries for A > ir/p are related to those for A < ir/p by reflection in 
the line A = ir/p. However, the phases themselves must be identified differently within the 
two regimes according to (f|). 

For D = oo the spins are restricted to discrete values 2imi/p where rii — 0,1, . . . p — 1 
and the Hamiltonian ([I]) becomes that of the p-state chiral clock model. The ground state is 
well known in this limit. For 0<A<7r/pitis ferromagnetic whereas for n/p < A < 27r/p, 
rij+i = Hi + 1, the increase in the chirality A favouring a twist in the spin ordering. At A = 
tt/p itself the ground state is infinitely degenerate, with any phase for which n i+ i — n,i = 
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or 1 for all i having equal energy. Such a point is often termed a multiphase point. It is 
expected on the basis of previous work that as D decreases from infinity the degeneracy 
will be lifted. Our aim is to explore the phase structure through an expansion in D -1 . We 
consider p > 3. 

To this end we require a notation capable of distinguishing the different phases stable 
at the multiphase point. Typically a stable ground state will consist of a sequence of bands 
where 7ij_i — = separated by walls with rij_i — rii = 1. [ £1^2 • • • 4n] will be used to 
describe a phase where the repeating sequence consists of m bands of length £1, £ 2 . . . l m . It 
may be helpful to list some examples for p = 6 

[12] ... I I 11 I 2 I 33 I .. . 

[23 2 ] ... I 00 I 111 I 222 I 33 I 444 | 555 | . . . 

[1] ... I I 1 I 2 I 3 I 4 I 5 I ... 

[00] ... 0000000 ... (5) 

where a vertical line is used to denote a wall. In the subsequent text we shall use the term 
/-band to describe a band of length I spins. For example [23 2 ] consists of a 2-band followed 
by two 3-bands. 

B. The Energy Differences 

The goal is to establish which of the infinite number of phases degenerate at the multi- 
phase point remain stable for finite D. This is done by using an expansion in inverse spin 
anisotropy, D^ 1 . The difficulty is that in order to check the stability of all commensurate 
phases certain terms must be calculated at all orders in D^ 1 . 

However, the relevant terms can be identified and the phase diagram constructed induc- 
tively using an argument first developed by Fisher and Selke to study the phase diagram of 
the axial nearest neighbour Ising model using a low temperature expansion^. We summarise 
their argument here. 
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Consider two phases [a] and [f3] which share a common boundary at a given order of a 
series expansion. Fisher and Selke showed that the first phase which can appear between 
them as the expansion is taken to higher orders is [7] = [a/3]. To check whether this phase 
does indeed appear the important energy difference is 

AE = 7l[ 7 ] Eft] - U[ a ] E[ a ] - 7l[/3] E[p] (6) 

where E[ a ] , E[p\ , E^ are the ground state energies per spin and ri[ a ], n\p\, n[ 7 ] the number of 
spins per period of [a], [f3] and [7] respectively 
There are three possibilities: 

(i) AE > and the boundary between [a] and \J3] remains stable to all orders. 

(ii) AE < and [a/3] appears as a stable phase in the vicinity of the [a] : [f3] boundary. 
The analysis must recommence about the new [a] : [a/3] and [a(3] : [/3] boundaries. 

(iii) AE = and [7] remains degenerate on the [a] : [/3] boundary to all orders. 

To explore the phase structure our goal is to calculate AE and investigate its sign. To 
this end we write 

Bi = 0- + Or (7) 

where 9® is the value of the spin 9i for D = 00 and expand the Hamiltonian ([!]) to quadratic 
order in the {9i}. 

n = n \ D=00 +^{M(0i-t - k + 4/4) 2 /2 + d%/2 - j(*?)7(2cf)} (8) 

i 

where 

cf = 006(^-6? + A), s^sinCd-^ + A). (9) 

We shall henceforth work with the quadratic approximation (||) to the chiral XY model. To 
leading order this gives the same results for the energy differences as the full Hamiltonian 

(I)- 

In equilibrium the energy of each phase must be minimal. Differentiating (Q) with respect 
to the 9i leads to the relation 
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D~6 t = J(c{ft_i - c% - clj, + 4 +1 B i+l + af - (10) 

which we shall need below. 

For the quadratic Hamiltonian 7i the energy differences can be calculated exactly. Let 
n [ a ] = n i and 7i[ 7 i = n and label the spins within the phases [a], [f3] and [7] by (ty, 7* 

respectively. Then, using (|8|) the energy of each phase relative to its value at D = 00 is 

m 

n x E [a] = E{Jcf (5<-i - o< + </cr) 2 /2 + #3?/2 - J(s?)7(2c?)}, 
i=i 

n 

(n - m)S M = E {Jcf(A_ 1 -A + S f/cf) 2 /2 + J DAf/2 - J( S f) 2 /(2cf)}, 

i=ni+l 
n 

n^ [7] = E{ JcJ^-i - 7 t + s7/c7) 2 /2 + £>7?/2 - J(,7) 2 /(2cJ)}. (11) 

8=1 

Because [7] = [a/5] we can choose to label the spins in such a way that 

7° = «°, l<i<n x , 

7° = A°, n a + 1 < i < n. (12) 

and take 

cf = cj, = 4, 1 < * < 

c f = c 7> s f = s 7> «i + 1 < i < n - ( 13 ) 

This means that we may drop the a, j3, 7 superscripts on the {q} and {sj} and the final 
terms in flllj) drop out when the energy difference is calculated 



AE = 

E[ J-MVi - 7i + Si/Q) 2 - Ci(5i_i - Qj + Sj/ci) 2 }/2 + L>(7 2 - fi 2 )/2] 



i=l 



+ E [ J{ci (7i_i - 7i + ^/q) 2 - q(A_i - A + Si/ci) 2 }/2 

i=ni+l 

+^(7 2 -A 2 )/2]. (14) 



This expression can be simplified considerably using (ID). Recalling the periodicity of the 
ground state phases which ensures 



&i = Oi ni+ i, (3i = (3 n - ni+i , 7j = 7 n+i V % (15) 

leads after some algebra to 

AE = Jci{(5 ni - /? n )(7i - 7m+i) - («i - Au+i)(7n - 7m)}/2. (16) 

Note that A.E depends only on the difference between a small number of spins. It is this 



which facilitates its calculation. The expression [16] is exact for the quadratic Hamiltonian |8| 
but only correct to leading order for the full Hamiltonian |TJ. However, this will be sufficient 
for the calculations presented below. 



C. Labelling the Spins 

The energy differences we are trying to calculate are independent of the labelling of the 
spins given that the conditions ( |T3|) hold. In general differences such as a ni — (3 n in the 
energy difference ( |T6"D will be polynomials in D~ l . Low order terms will cancel when the 
difference in ( pf ) is taken, in such a way that the final result becomes independent of the 
labelling. However calculationally the problem is simplified by a careful choice of spin labels 
which allow the leading order contribution to AE to be obtained directly. 

We first point out that every commensurate ground state of the Hamiltonian (1) has 
two points of mirror symmetry in each period evenly spaced along the chain!. It is possible 
to distinguish two cases. For states of odd period half the symmetry points are located 
on lattice sites and the other half between lattice sites. By symmetry the spins located on 
the lattice sites corresponding to mirror symmetry points do not deviate from their D = oo 
position for finite D (6 = 0). For states of even period the symmetry points are located either 
all between or all on lattice sites. For states formed by the branching process [a] + \J3] ^> [af3] 
only the former can occur. 

As the branching process [a] + [f3] =3- [7] = [a(3] proceeds states are made up in two ways 

(i) odd + odd even, 

(ii) odd + even =^> odd. 
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(A moment's reflection shows that even + even =>- even never occurs because, as the ground 
states are formed inductively, no neighbouring even phases ever appear.) 

Our aim is to choose the origin for the labelling of the spins in such a way that the 
formula flip] ) for the energy difference is simplified. It is necessary to consider each of the 
two cases (i) and (ii) separately, 
(i) odd + odd =>- even 

We recall the notation [a] + [/3] =>- [7] with [a] and [7] having n\ and n spins respectively in 
each period. A convenient choice of labelling is 

< =0; P n = 0. (17) 

It will be useful later to write the spin deviations within a period of each phase in a way 
that explicitly displays the symmetry 

{O-i} = Oi 2 ■ ■ . 5( m -l)/2, <S(ni+l)/2 • • • <5ni-2, "m-l, «m } 

= {&i, a 2 ■ ■ ■ 5(ni-l)/2, — «(ni-l)/2 • • • — «2, — &l, 0} (18) 

{Pi} = {Prn+1, Pm+2 ■ ■ ■ /3(n+ni-l)/2, +l)/2 • • • Al-2, /?n-l, 

— {Prn+1, Pm+2 ■ ■ ■ P(n+ ni -l)/2), ~ P{n+ ni -l)/2, ■ ■ ■ ~ Pm+2, ~ 0}. (19) 

(a) and (/3) combine to give an even state (7) with symmetry points between i = {n\ — l)/2, 
% = [m + 1) /2 and % = (n + n x - l)/2, i = (n + n x + l)/2. 

{7i} = {71,72 ■■■ 7(ni-l)/2,7(m+l)/2 ••• 7n 1 -l,7n 1 ,7n 1 +l,7ni+2- • ■ 
7(n+m-l)/2,7(n+ni+l)/2 • • • 7n-l> 7n} 
= {7l, 72 ■ ■ ■ 7(ni-l)/2, -7(m-l)/2, ■ ■ • ~ 72, -Jn, ~7n-l • • • 

— 7(n+ni+l)/2> 7(n+ni+l)/2 • • • 7n-l, 7n}- (20) 

Because a ni — (3 n = 0, ( ]TB| ) simplifies immediately to 
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AE = -Jcxi&x - /3 ni+ i)(7n - TnJ/2. (21) 

(ii) odd + even =>• odd 

We take [a] odd and [f3] even. Choosing ar ni+ x)/2 = 

[a i } = {a ll a 2 ■ ■ ■ fi( ni _i)/ 2 , 0, -a( ni _i)/ 2 , • • • - a 2 , —ax}. (22) 

A consistent choice of labelling for [j3] which results in the correct final state is to take the 
mirror symmetry points to lie between spins i = and i = 1 and between i = (n — n{)/2 
and i = (n — n 1 )/2 + 1. Hence we may write 

{Pi} = {Pni+li An+2 • • • An+m)/2, — P(n+m)/2 ■ ■ ■ — Pni+2, ~ Pni+l}- (23) 

It is immediately apparent from (|22|) and ( p3|) that 

am ~ Pn = ~(ax - P ni+1 ). (24) 

Using [a] and to construct [7] will preserve a point of mirror symmetry at i = (ni + l)/2. 
Hence 

{7i} = {7l, 72 ■ ■ ■ 7(ni-l)/2, 0, -7(m-l)/2 • • • - 72, -7l, -7n, -7n-l • • • 

-7(n+ni)/2, 7(n+m)/2, • • • 7n-l, In} (25) 

from which it follows that 

7n - 7m = -(7m+i - 7i)- (26) 
Using (|24]) and (p6|) the energy difference (^) simplifies to 

A£ = -Jci(a! - /9 ni +i)(7n - 7m)- (27) 



D. Recursion equations 



The next step is the calculation of the spin deviations in the formulae fl21"D and (27|). To 
do this we start from equation (^) which followed from minimising the ground state energy. 
Let 
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(28) 



Substituting into (10) and equating like powers of D 1 gives 

9\ = J( Si - s i+l ), (29) 

C = J(qC/ - - Ci+i^ -1 + Q+i^+i 1 ), n>\. (30) 

Again it is necessary to consider separately the combination of states with different symme- 
tries. 

(i) odd + odd — » even 

Note firstly that a] — Pn i+i depends only on the {s"} and {sf }, that is only on the value of 
the spins for D = oo. Let 

«i - #t 1+ < = 0, i < no, (31) 
a- - = ao, i = no. (32) 



Then a consequence of the spin labelling and symmetry summarised by equations (18) and 



(|T|) is that 

a^- 01^ = 0, t<n , (33) 
atli - Pn x -i = -ao, i = no- (34) 

It is apparent from the recursion equations ([3(]) that after one step of the iteration the 
spin differences cij — $ ni +i and — $m-i with i = (no — 1) will attain a non-zero value; after 
a second step these differences with % = n — 2 will become non-zero and so on. Therefore 
iterating (n — 1) times gives the leading order result 

&i - K+i = J 110 ' 1 cosiir/vY^av/D^. (35) 

Note that we have assumed that q = cos(7r/p) V %. This is true to leading order near 
A = n/p. 
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It is also appropriate to mention here that further iteration of the linear equations will 
generate corrections to fl35|) which are correct for the quadratic Hamiltonian (||) but not for 
the full chiral clock model ([!]). However, these terms are not necessary for our argument. 

As a consequence of the conditions (|1^) and because [7] = [a [3] it must also hold that 

7; - = °> ~n <i< n , (36) 
7*-7n 1+ i = ao, i = ±n . (37) 

Iterating n times using ([$(]) gives 

% - 7 m = 2 J n ° cos(vr/ p) n °a /D no+1 (38) 

where the factor 2 appears because terms which iterate both from the right and left along 
the chain contribute to this order. Hence the energy difference (^1]) is 

AE = -J 2no cos(n/p) 2no al/D 2no+1 . (39) 

(ii) odd + even =>• odd 

We assume as before that (|31|) and ([32]) hold. It follows immediately from considering (|22|) 



and (E3) that 



<* t - - Pn 1+ i = 0, -(n - 1) < i < n , (40) 

"no _ Pni+no = °0> a -(n -l) ~ Ati-(no-l) = ~ a °- ( 41 ) 



Similarly 



7i - 7n 1+i = 0, -(n - 1) < i < no, (42) 
ll ~ lm+i = «o, i = n ,-(n - 1). (43) 

Iterating (uq — 1) times 

fix - = J"'- 1 cos(7r/p) ?1( - 1 a / J D" (44) 

from terms iterating from the right and 

12 



7n - 7m = J* 10 " 1 cos(7r/p)"°- 1 ao/ J D no (45) 
from terms iterating from the left. Using (|27D the energy difference is 

AE = -J 2 «o-i cos{-K/p) 2na - l al/D 2no . (46) 



Both of the energy differences fl39|) and (f46j) are negative. Hence for p > 3 all phases 



appear for large .D near the multiphase point in agreement with the conclusions of Chou 
and Griffith!. 

Noting that Sj = ± sin(7r/p) V i it is apparent by inspection that a 2 , = 4J 2 sin 2 (7r/p). We 
prove in Appendix A that = [(nr 7 i + l)/2] — 1 where [n] is the integer part of n. Therefore 
for all combinations of phases 

AE = - J n H" 2 cos(7r/p) ri M- 2 a 2 / J D n H- 1 . (47) 

Correction terms 0(1/ D n w) will arise from contributions to the energy difference (|T6|) and 
the recursion equations (|T0D from non-harmonic terms in the Hamiltonian dH). There will 
also be correction terms within the harmonic approximation (H) itself arising from further 



iteration of the recursion equations (|T0D. 



E. Phase Widths 

An advantage of the formalism presented above is that it allows a calculation of the 
widths of the long-period phases. We define A a7 , A^, and A a/3 by 

E [a] (A a7 ) =E [7] (A trr ), (48) 
E m (A^) = E b] (A^), (49) 

E[ a )(A a p) — E\p](Aap). (50) 

For a stable phase [7] A a p will lie between A ai and A^ 7 and if the phase [7] is not too wide 
a Taylor expansion gives 
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E[ a ](A ai ) = E[ a ](A a p) + E'^(A ai — A a p), (51) 
E[ 7 ](A a j) = Efa](A a p) + E^(A ay — A a p) (52) 

where ' denotes a derivative with repect to A. 



Subtracting ( p2|) from (|5l|) and using (|4q) gives 



(^[a] _ -£[ 7 ])(Aa7 - A Qj9 ) = ^ [7 ](A Q/3 ) - E[ Q ](A a/3 ) (53) 

= £ [7] ( A aP ) - ^£ [q] ( A aP ) - (A a/3 ) (54) 

n [7] n [7] 

- ^ (55) 

7l [7 ] 

where in the penultimate step we have used ( (5"0| ) and in the final step the definition of AE, 
equation @. Writing down a similar expression for A^ 7 — A a @ and combining it with ( 55|) 
gives an expression for the width of the phase [7] 

W b] = (A^ - A Q7 ) = AE{(E' m - E'^)- 1 - (E' [a] - E' {l] )- l }/n b] . (56) 

For D = 00 the energy per spin of a ground state phase at the multiphase point [a] say 

E [a] = - J{l [a ] cos A + (1 - l [a] ) cos(-2tt/p + A)} (57) 

where l[ a ] is the fraction of nearest-neighbour ferromagnetic bonds. Differentiating the ex- 
pression (l57f) and similar formulae for [0\ and [7] and substituting into (|56f) gives to leading 



is 



order 

^ = Aff(/ M - jj^) 

It1 2Jn H sin(7r/p)(/ [/3 ] - Z[ 7 ])(Z[ a ] - Z[ 7 ]) ' 

It is not hard to show inductively that 

l[/3] ~ l[a] = (ri[ a ]n[/3]) -1 , Z[/3] - Z[ 7 ] = (n^Tlf,])" 1 , l[ a ] ~ i[ 7 ] = (n [a ]?2[ 7 ]) _1 . (59) 

Therefore using fl47D we finally obtain 

W M = -2J n H- 1 cos(7r/p) n H- 2 sin(7r/p)n [7] / J D n w- 1 . (60) 
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III. XY MODEL WITH COMPETING FIRST-AND SECOND-NEIGHBOUR 
INTERACTIONS AND P-FOLD SPIN ANISOTROPY 



A. Definitions and Notation 

In the second part of the paper we extend the formalism developed earlier to obtain new 
results for a more complex spin model. This is the XY model with competing first- and 
second-neighbours interactions and p-fold spin anisotropy. Each classical XY spin vector 
lies in a plane perpendicular to the z-axis and has unit length. The Hamiltonian can be 
written 

H = £ \—J\ cos(^_! - Oi) + J 2 cos(^_ 2 - 9i) - D(cos(p0i) - l)/p 2 } (61) 

i 

where $i is the angle between the spin located at site i and, say, the x-axis. Competition is 
introduced along the ^-direction by taking the first and second neighbour interactions to be 
ferromagnetic and anti-ferromagnetic, respectively {J\ > 0, J 2 > 0). x = J 2 /J\ will prove 
an important variable in the description of the phase diagram. 

The parameter D > models a p-fold spin anisotropy in the (x, y) plane. The ground 
state in the two limits D = and D = oo is well understood. For D = it is ferromagnetic 
for x < 1/4. For x > 1/4 it exhibits helical order with a wave vector q = qz which is, in 
general, incommensurate with the underlying lattice. The magnitude of the wave vector is 
determined by the exchange energies through the relation cosg = (4x) _1 . 

For D = oo, however, the spin angles 6i are constrained to take one of the discrete set of 
values 2iini/p, rii — 0, 1, . . .p—1. The Hamiltonian (|6l|) then reduces to ap-state clock model 
with competing interactions. The ground state now has a very different character: only a 
few short-period commensurate phases are stable as x is varied. Boundaries between the 
different ground states can either be simple first-order transitions with only the neighbouring 
phases being stable or multiphase points at which an infinite number of phases have the same 
energy. 
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To fully describe the ground states and the phases degenerate at the multiphase points 
it is necessary to extend the notation introduced in Section [11 A[ Note that at multiphase 



points such as x = 1 for p — 6 all phases with rii — rii-i = 1, 2, with the proviso that 
rii + i — rii = rii — rij_i = 2 is not allowed, are stable. Hence the natural definition of a wall is 
as lying between spins i and % — 1 for which — n,i_i = 2. The term band is used as before 
to describe the sequence of spins between two walls. The phases stable at the multiphase 
point can then be more easily described as those containing only bands of lengths > 2. 

A given state will be labelled by < £\, £2 ■ ■ ■ £ m > where the repeating sequence comprises 
bands of length £\, £2 . . . £ m . It may be helpful to list some examples for p = 6 

< 2 > ... I 01 I 34 I 01 I .. . 

< 23 > ... I 01 I 345 I 12 I 450 I .. . 

< 3 2 4 > ... I 012 I 450 I 2345 | 123 | 501 | 3450 | . . . 

<oo> ...01234... (62) 

Compare the notation using square brackets [€i,£ 2 . . .£ m ] introduced in Section [11 A| where 



the walls correspond to — n^-i = 1 in a background matrix of rii — rij_i = 0. This will 
also be needed here. 

We are now in a position to describe the ground states of the Hamiltonian ( |6lD for all 



values of p and D = 00. The results which were obtained by comparing the energies of the 
possible ground states and checked using the Floria-Griffiths algorithm! to ensure no states 
were missed are summarised in Figure 1. 

(i) For p = 2 and 3 [00] is stable for x < 1/2 and [2] for x > 1/2. x = l/2isa multiphase 
point. 

(ii) For p = 4 [00] is stable for x < 1/2 and < 00 >= [1] (together with the phase . . . 002200 
. . . ) for x > 1/2. Again x = 1/2 is a multiphase point. 

(iii) For p = 5 for x < xf^ = {1 + cos(27r/5)} _1 /2 [00] is stable. For x > xo, < 00 > is 
stable, xq is not a multiphase point. 

(iv) For p > 6 there is a common trend for small x. For x < Xi = {1 + cos(27r/p)} _1 /2 
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[oo] is stable. For x-f < x < Xq = {cos(27r/p) — cos(67r/p)}/2{cos(47r/p) — cos(67r/p)} the 
phase < oo > appears. For x > Xq < 2 > is stable. (For p = 6 the state . . . 003300 ... is 
degenerate with < 2 >.) For p = 6 and 7 the phases listed above provide all the ground 
states. For p > 8 other phases are stable for higher values of x, but these phases will not 
concern us here. At the points the phases [oo] and < oo > coexist but no other phases 
are stable and there is a first order transition. Xq however is a multiphase point where all 
phases comprising bands of length > 2 are degenerate. 

Our aim is to apply an expansion in D^ 1 to understand the phase structure in the vicinity 
of the multiphase points xffl . We proceed in the same way as in Section |T|. The results are 
listed in Section [III F| . A brief account of this work has appeared elsewhereB. 



B. The Energy Differences 



Our approach follows that described in Section [II B| for the chiral clock model. The first 
step as before is to expand the Hamiltonian flSl| ) to second order in the 9i, the deviations of 
the spins from their positions at D = oo, which were defined by equation (|7|). This gives 

H = H |d=oo +XX*^,[i](^-i + A i,[i]) 2 ~~ M s l[i}) 2 / C i[i] 

i 

+ ^« [2] )7<[ 2 ] - - ~9i + Aj [2] ) 2 + D§n/2 (63) 

where 

S i,[m] = S i n (^i°-[m] — ^i)) C i,[m] = COS (@i-[rn} ~ 

A» _ J I J jd _ J J> (RA\ 

^-\[m] — °i,[m]l H,[m\i u i,[m] ~ °rn^i^ m y \v<±) 

Minimising (|63|) with respect to 6i leads to the equation 

D9i = Ji,[i](0t-i -0i + - Ji+i,[i]{0i - 9 i+ i + A j+1)[1] ) 

—Ji,[2]{0i-2 - 9i + A i>[2 ]) + J i+ 2,[2]{9i - 9 i+2 + A i+2)[2] ). (65) 

It is possible to calculate the energy difference AE, defined by (|5]), exactly for the 
quadratic Hamiltonian (|63|) . As before we consider the branching process [a] + \J3] =^> [a(5] = 
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[7], let n[ a ] = ni and n[ 7 ] = n and identify the spins within the phases [a], [/3] and [7] by 
Qi, Pi and 7j respectively We choose to label the spins such that 

7° = l<i< m, 

7? = Pi m + 1 < i < n. (66) 

from which it follows that 



A" 

i,[m] 


— A 7 


TO — 7" 7 
"^[m] ~~ i,{m) 


A 13 


- A 7 


jP _ t7 
J i,\m] ~ J i,\m) 



The energy difference (Bf) may now be written 



(67) 



2AE = X){[J»,[i](7i-i - 7i + ^,[i]) 2 - ^,[i](«i-i - «i + Ai i[2 ]) 2 ] 

- [Ji,[2](li-2 - 7i + A i,[2]) 2 - ^,[2](«i-2 - «i + A i,[2]) 2 ] 

+ £(7?-« 4 2 )} 

n 

+ E {[Ji,[i](7i-i - li + A 4 ,[i]) 2 - - A + A iiW ) 2 ] 

- [Ji,[2](li-2 ~ li + ^i,[2]) 2 ~ Ji,[2](Pi-2 ~ Pi + Aj 5 [2]) 2 ] 

+ D(j!-PI)}- (68) 



Using equation (|65|) this can be simplified to give 

2AE = Ji,p.]{(a ni - /9 n )(7i - Tm+i) ~ («1 _ Au+i)(7n - 7m)} 

- Ji,[ 2 ]{(a ni _i - /3n-i)(7i - 7m+i) - («l ~ /Wi)(7n-i - 7m-l)} 

- ^2,[2]{(a m - /3„)(72 - 7n 1+ 2) - («2 ~ /3ni+2)(7n - 7m)}- (69) 

Again the energy difference depends only on the difference of the spin deviations on a few 
sites. Given a careful labelling of the states these can be calculated to leading order and 
their signs determined. 



18 



C. Labelling the Spins 

We label the spins in a way identical to that described in Section |11 Q This simplifies 
the formula for the energy difference AE as follows: 
(i) odd + odd =>- even 

With the choice ( |17|) it can be read off immediately from (18), (|19|) and ( pOl) that 



C^ni Ai 0) 

OCl — Au+1 = ~ (<5ni-l — Ai-l)> 

(7i - Tm+i) = (7« - 7m) (70) 



and hence from (|6 



2A£ = - Ji,[i](5i - A»i+i)(7n - 7m) 

+ ^i,[2](ai - ^ m+ i)(7n-i - 7m -1 + 7i - 7m+l) 

+ ^2,[2](«2 - A»i+2)(7n - 7nJ- (71) 



fii) odd + even =>- odd 



Using (H), (H), and © 

— Ai-1 — ~ (^2 — Aii +2)) 
«ni - /3 n = -(«i - A tl+ l), 
(7„_1 - 7m-l) = -(7ni+2 - 72), 

{in ~ 7m) = -(7ni+l - 7l) ■ (72) 



Therefore 



2AE = 2J 1)[1] (5 1 - /3 ni+1 )(7! - 7 m+1 ) 

+ «7l,[2]{(«2 - Au+2)(7l - 7ni+l) + («i - An+l)(72 - 7m+2)} 
+ ^2,[2]{(«1 - An+l)(72 ~ Tm+a) + (& 2 - A11+2X71 - Tm+l)}- 

(73) 
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D. Recursion equations 



We follow Section [11 D| in using recursion equations to derive the leading order terms 
in the energy differences ( |7T|) and (^). Substituting fl28|) into the energy minimisation 
equations flSo] ) gives 

9} = Ji(s i} [i] - Sj+iji]) - J 2 (Sj,[2] - S i+2) [2]) (74) 
~ Ji,[2](Qi-2 ~ &i~ ) + Ji+2,{2}(0?~ 1 - (75) 

For the model with second-neighbour interactions it is necessary to consider four different 
cases when calculating the energy differences, 
(i) odd + odd =>- even; even starting position 
Let 

a 2n ~ Pni+2n = a 0; a -2n ~ A^-2n = _a 0, (76) 

that is i is even when a\ — +i is first non-zero. We consider uq > 1 throughout. This 
implies 

7i - 7ni+i = 0, - 2n o < * < 2n , 

7i - 7ni+i = «o, i = ±2n . (77) 

The recursion equations ([7^) show that the spin deviations (5j — I3 ni+i for i = ±(2n — f) 
and i = ±(2n -2) will be C(f/£> 2 ); for i = ±(2n -3) and i = ±(2n -4) will be 0(1/ D 3 ) 
and so forth. Hence 

& x - P m+1 ~ 0(f/.D no+1 ), 7n - 7m ~ O(l/ J D n0+1 ), 
^-% 1+1 ^O(l/D n0+1 ), a 2 -P ni+2 ~0(l/D n °) (78) 

and the leading order term in the energy difference (ffll) is 

20 



2AE = J 2)[2] (a 2 - P ni+2 )(jn ~ 7m). (79) 
From the recursion equations ( |75|) iterating (no — 1) times 

a 2 - K+2 = {a {-J 2 ) n °-\lf no ^cir l ~ nw[2n ° ,2] }/D no (80) 
where n w (i,j) is the number of walls between i and j and 

cos(a;°_ 2 — a°) = c 2c = cos(67r/p) (81) 
if the iteration step jumps across a wall and 

cos(o;°_ 2 — a°) = c 2 j = cos(47r/p) (82) 

if it does not. 

Similarly iterating n times and collecting equal terms from the left and right 

7n - 7 m = {2a (-J 2 ) no c 26 42n0i0] c 2 r n42rio ' 0] }/ J D no+1 . (83) 
Substituting (|(]) and (|83l) into (^|) gives 

_ {J 2[2]a 2 (-J 2 ) 2 ^-l c ^[ 2 ^' 2 ]+™»[ 2 ™0,0] c 2no-n IU [2no,2]-n tl ,[2r i .o,0]-l^ D 2no+l 

a n = 0. Therefore the band containing i = n of the state < a > must be of odd length. 1- 
bands are forbidden. If the band is of length 3, J 2 ,[2] = <^2C 2c and n w [2n , 0] = n w [2n , 2] + 1. 
If it is of length > 5, J 2) [ 2 ] = J 2 c 2 j and n w [2n , 0] = n w [2n , 2]. In both cases 

AE = _{ a 2 (J 2 ) 2 «o c 2n»[2no,0] c 2n -2n l „[2no,0] }/jD2no+ l (g5) 

where the relevance of the subscript on AE will become apparent later. 

(ii) odd + odd =>- even; odd starting position 

Let 

a]-(3l i+l = 0, -(2n - 1) < i < (2n - 1), 

a 2n -l ~ Pni+2n -l = a 0> a -(2n -l) — Pn 1 -(2n Q -l) = ~ a (86) 
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where we now consider i odd for a] — f3\ 1+i first non-zero. This implies 

7--7^ 1+i = 0, -(2no-l) <i<(2no-l), 
li ~ ln 1+i = °o, i = ±( 2n o - 1). 

Iterating the recursion equations ( |75| ) shows that 



no+l\ 



at - (3 ni+1 ~ Oil/I?*), % - 7n: ~ 
7i - 7n 1+ i ~ O(l/£> no ), « 2 - ^ ni+2 ~ 0(1/L> n °) 

and hence the leading order term in the energy difference (|7l| ) is 

2A£ = Ji,[ 2 ](ai - /3m+i)(7i - 7m+i + 7»-i ~ Tm-i)- 
From the recursion equations (|75|) iterating (no — 1) times 

«1 - Aii+1 = 7l - 7m+l = 7n-l ~ Jm-l = 

r„ i j \no— i „ n K>[ 2n o— no— i— »*T»[2no— l,i] i / n n o 
•K\2] = ^2C 2 i as the band containing the spin i = n must be of length > 3. Therefore 

AE 2 = { a 2 (j2) 2no-l c 2n ra [2no-l,l] c 2n -l-2n 1 „[2no-l,l]| /jD 2no > 

(iii) odd + even =>- odd; even starting position 
Let 

«■ - = 0, -(2n - 1) < % < 2n , 

a ln ~ Pni+2n = fl 0> a_(2n -l) — Aii-(2na-l) = ~~ a 0; 
2n +l ~~ 0m+2n o +l = a l> a -2n ~ Pm-2n = — fl l 



ft 



and similarly for the 7* — 7^ 1+i . Noting from the recursion equations (|75|) that 

fix - ~ 0(1//J" O+1 ), 7i " 7m+i ~ 0(1/D n ° +1 ), 

72 - 7n 1+2 ~ 0(1/D^), a 2 - fi ni+2 ~ 0(1/1}"°) 
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and that whatever the arrangement of walls symmetry implies J\m = J2J2] the energy 
difference (F3D becomes 



A£ = Ji,[ 2 ]{(5 2 - P ni +2)(jl - Tm+l) + ("i - Au+l)(72 - Tni+2)}- (94) 

Iterating the recursion equations (JF|) (no — I) times 

~ n ~ ~ ft j \n -l n w [2n ,2] n -l-n w [2n ,2]^ , nno / nt .s 

"2 - Pm+2 = 72 - 7nj+2 = {a (-J 2 J c 2c c 2i \/D °. (95) 

Calculation of «i — /?„ 1+ i is slightly more involved as iterating n times contributions can 
arise either from hops from i = 2n which include a Ji term or from hops from i = 2n$ + 1 
which include only J 2 terms 

a r / t \na— 1 t n m [2no , 1] — 1 no-n w [2no,l] — 1 / / ~ \ i ~ \ 

oil ~ Pm+i = |a (-^2) J\c 2c C2i {C2cCu{n - n ) + c 2 iCi c n ) 

i „ l T \n n w [2n +l,l] n -n w [2n +l,l] 

+ <M — J 2> c 2c c 2i 

- ^(-Ja)" ^; 1 "^ - 1 ^ 11 ^"" 4 "^ 0-1 ^ 11 }/^^ 1 (96) 

where n is the number of paths in which the J\ hop crosses a wall between 2i%q and 1, 
Cic = cos(47r/p), and Cu = cos(27r/p). The expression for ji — 7„ 1+ i is similar but the last 
term, as it corresponds to iteration from the left, contributes with opposite sign. Therefore 
the energy difference (|73|) is 

AE 3 = 

{2J m al(J 2 f^-\j lC ^ - n ) + c 2l c lc n ) 

OT „ „ / T \2nn-l ™™[2no,2]+n m [2no+l,l] 2no-l-n„,[2no,2]-n„,[27io+l,lh / n2n.n+l 
— ^l,[2] a a l \ J 2) C 2c c 2i } I U 

(97) 

(iv) odd + even =^> odd; odd starting position 
Let 

a] - f3l i+i = 0, -(2no - 2) < i < 2n - 1, 

a 2n -l — /^ni+2n -l = a 0; a -(2n -2) — /^ni-(2n -2) = ~ a ; 

a 2n _ /^ni+2n = °1> a -(2n -l) — /^ni-(2n -l) = _ fl l (9§) 
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and similarly for the 7* — 7 ni+i . Iterating the equations (|75|) (no — 1) times 

«1 - Pn 1+ 1 =11- 7n 1+ l = W- J 2 )«0-l c ^[2no-l,l] c no-l-n»[2no-l,l] }/jD nO ; (gg) 

"2 - /^ni+2 = 72 - 7ni+2 

r / r Min-2 t n,„[2n — 1,21-1 no— n«,[2riQ— 1,2]— 2 / / ^ /\ / \ 

= |a (-J 2 ) ^ic 2c c 2 i (c2 C Cii(r2 - 1 - n ) + c 2i Ci c n ) 

+ai(-J 2 ) no ~ 1 C2r [2no,2] c 2 t r 1 " nt " [2no ' 2] }/^ no (100) 

where n' is the number of paths in which the J\ hop crosses a wall between 2no — 1 and 2. 
Using these formulae the energy difference (|73|) is 

AEi = 

r t 2/ t \2nn-2 2n„[2n -l,l] 2n -2-2n u ,[2n -l,l] 

{— Ji,[i\a {J2) c 2c c 2i 

o t 2/7 \2nn-3 7 n w [2no-l,2]+n TO [2no-l,l]-l 2no— 3— nto[2no— 1,2]— nro[2no— 1,1] 
—ZJl,[2]a [J2) JlC 2c c 2i 

(c2cCu(n - 1 - n ) + c 2i Ci c no) 

1 07 ^ ^ f 7 N2n -2 n w [2n -l,l]+n w [2n ,2] 2n -2-n u ,[2n -l,l]-n m [2n ,2] -1 / n 2nn 

+zj l [ 2 ]a a 1 [J 2 ) c 2c c 2i )ju 

(101) 



E. Phases bounding (00) 

The formulae for the energy differences calculated in the previous section may no longer 
hold when one of the initial phases is (00) = . . . 0123 . . .. Therefore this case is now treated 
separately As the phase diagram is constructed recursively the phases which will appear on 
the (00) boundary are (m), m — 3, 4, 5 . . .. Hence the energy differences we wish to calculate 
are 

A£<°°> = (m + l)£ (w+ i> - mE {m) - E (oo) . (102) 
In the (00) phase the spins remain in their clock positions as D is reduced from 00. 
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We must again consider four possibilities. 

(i) (An + 3) + (oo) =» (An + 4) 

We consider n > 1 throughout. Choose (f3) = (oo). Then (3 = and n\ — n — 1. The initial 
conditions are 

aj = 0, -2n < i < 2n, (103) 

«L = ~a l -2n = a o, (104) 
7^ = 0, -(2ra + l) <i<2n, (105) 

72n = -7i(2«+i) = fl o- (106) 



Noting that j ni = the calculation follows that described in Section p[~D|(i). The energy 
difference is given by equation ([55]) with n w [2n, 0] = 

AEt ] = -{al(J 2 ) 2n c%}/D 2n+1 . (107) 

(ii) (An + 1) + (oo) g> (4ra + 2) 

Again choose (/3) = (oo). The initial conditions are now 

a} = 0, -(2n- 1) < i < In- 1, (108) 
= -«-(2n-l) = fl o, (109) 

7^ = 0, -2n<i<2n-l, (110) 

7an-l = -l X -2n = a 0- (HI) 

Recalling that n\ —n — 1, 

7i ~ lm-i ~ «i ~ a 2 ~ C(l/-D n_1 ), 

7n-i ~ J ni -In- ~ (112) 

and that the symmetry of (7) implies 

7i = -7m-i (H3) 
the leading term in the energy difference ([H]) is 
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AE = 2 Ji,[ 2 ]fii7i- (114) 
Iterating the recursion equations ( f75| ) (n — 1) times to obtain and 71 gives a result 

A£ 2 <00> = {a^^) 2 ™- 1 ^- 1 }//} 2 ™. (115) 

(iii) (An + 2) + (00) =» (4n + 3) 

To use the energy difference formula (60) ((3) must be even. Therefore we choose (a) = (00) 
implying ni = 1 and a, = V z. The initial conditions are 

/£ 1+i = 0, -(2n-l)<i<2n, 

/^ni+2n = — Aii-(2n-l) = fl An+2n+l = ~Pni-2n = fl l 

7-=0, -(2n-l) <i< (2n+l), 

7-L+l = ~7i(2n-l) = a 0. 7<L + 2 = ~l-2n = °1" (H 6 ) 

Using the symmetry properties of ((3) and (7) 

/3 ni+i = -(3 ni - i+ i, 7 = -7- i+2 (117) 

and noting that 

h ~ 73 ~ ©(l/i^ 1 ), & ~ 72 ~ 0(1/^0") (118) 
the leading order contribution to the energy difference (|73|) is 

2AE = (J 1>[2] + J 2 , [2 ])(/3 3 72 + fa%). (119) 
Iterating the recursion equations (|7^) to obtain the spin deviations the final result is 

AEt ] = Jl n - l c 2 £- l {al{2J x nc u + J 2 c 2i ) - 2a a 1 J 2 c 2i }/D 2n+1 . (120) 

(iv) (An) + (00) =» (An + 1) 

Choosing (a) = (00) the initial conditions are 
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/£ 1+i = 0, -(2n-2)<i<2n-l, 
7 i=0, -(2n - 2) < i < 2n, 

lL = -^-(2n~2) = a 0, lln+l = -J-(2n-l) = °1- ( 121 ) 



Using the symmetry properties (117) the leading order contribution to the energy difference 
(0) ^ 

2A£ = 2J li[1] &7 2 + (J 1>[2] + J 2j[2 ])(^372 + /3 2 73 - £272) • (122) 
Using the recursion equations (|75|) we obtain 



A£f°> = - J 2 2 "- 2 C 2r 2 {a2(j lCli (2n - 1) + J 2 c 2j ) - 2a a x J 2 c 2i }/D 2n . (123) 

There will be correction terms to the formulae for the energy differences which arise both 
from non-harmonic terms in the Hamiltonian floTD and from further iteration of the recursion 
equation Q75p. These will be more dangerous than for the chiral XY model as they may carry 
an additional factor n^/D where the n 2 , comes from say dividing a J 2 step into two J\ steps 
and placing them at any position along the chain. Therefore for large no these terms could 
dominate, possibly changing the sign of the AE and the sequence of phases could terminate 
for any finite D. Analytic calculation of the correction terms would be prohibitively difficult. 
However numerical results show no sign of deviations from the leading behaviour for nr 7 i=15. 
A similar mechanism has been desribed for the ANNNI model at finite temperatures^. 

To apply the formulae for the energy differences it is necessary to ascertain which AE^ 
and AE^°\ k — 1, 2, 3, 4 should be used at each step of the iteration procedure. By checking 
low order examples or by an inductive argument similar to that given in Appendix A one 
finds that AE^ and AE^ are relevant for a final phase with nui = k [mod4]. 

The values of a and ai 

a = J2{sin(67r/p) — sin(47r/p)}/D 

a-y = -J^sin^Tr/p) - sm(2n/p)}/D (124) 
follow immediately from (|75|). 
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F. The Sequence of Phases 



The stability of the short-period phases is most easily checked by an explicit evaluation 
of their energy using the Hamiltonian fl63|). Phases which appear 0(1/D) and 0(1/D 2 ) are 
shown schematically in Figure 2. Stable boundaries are denoted by a vertical line; boundaries 
at which an infinite number of phases still coexist 0(1/ D 2 ) by a star. By using the energy 
differences derived above it is possible to calculate which of these phases are in fact stable 
when higher order terms in 1/D are considered. The results depend sensitively on p and so 
we consider each value in turn. 
p = 6 

For p = 6 the calculation of terms 0(1/ D 2 ) establishes that the (2) : (3) boundary is stable. 
However all other boundaries are multiphase lines and therefore any phase which contains 
only bands of lengths > 3 may appear in the phase diagram. To understand which of these 
phases are stable we need the signs of the energy differences. 

The energy differences are simplified by noting that for p — 6, ai = 0. To obtain the sign 
of AEi given by equation ( |101|) two cases must be considered. If there is a wall on the axis 



of symmetry (i.e. between «o and a{) then, recalling that all bands are of length at least 
3, Jim = J\Ci C) Ji,[2] = J2C2C, and n w [2n — 1,2] = n w [2n — 1, 1]. Putting in values for c 2c , 
etc. gives 

AEi = 3J 1 J 2no (l/2) {2n °- 2n ™- 1} (5/4: - n + 3n' /2)/D 2n °. (125) 

If there is no wall on the axis of symmetry J^pj = J\Cu, J^p] = Ji^n-, and n w [2n — 1,2] = 
n w [2n — 1, 1]. Therefore 

A E ± = 3J 1 J 2no (l/2) {2no ~ 2n - +1} (l - 2n + 3n' )/D 2n °. (126) 

which is negative. Using similar arguments, or by inspection, all the other energy differences 
are found to be negative0. Therefore all phases which contain only bands of lengths > 3 
appear in the phase diagram. 
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p = 7 

For p = 7 it can be shown, using arguments similar to those presented for p = 6 that all the 
energy differences are negative for large enough n . The cases where n is small and n' Q and 
no are close in value to no must be checked independently. One finds that all the relevant 
cases are negative except AE\ for no = 2, n' — 1. This determines the stability of (2 3 3). 
Therfore for p = 7 the (2) : (223) boundary is not split. All phases lying between (223) and 
(oo) are however stable. 
p = 8 



The energy differences AE , AEi, AE 2 , and AE 3 given by equations (|85|), (|101|) , (|9~1~D , and 



(|9~7|) are zero. Therefore higher order terms which are prohibitively difficult to calculate are 
needed to establish the signs of the energy differences. 

However, it is possible to obtain an expansion of the energies of the low order phases 
numerically directly from the Hamiltonian fl61~l). (The quadratic approximation (|63|) to the 
Hamiltonian may not be sufficient to pick up the correct leading behaviour in this case. 
Comparing these energies one finds that at least all phases expected to appear 0(1/ D 5 ) 
(n( 7 ) < 12) are stable. 
p = 9 

For p = 9 the situation is complicated and no clear pattern of phases emerges. Results 
0(1/ D 2 ) indicate that only phases which are made up of bands of lengths < 4 can appear. 
One finds by inspection AE < (equation (|85|)); AE 2 > (equation (0)). For large enough 
n if there is a wall on the axis of symmetry AE 3 < and AE\ > from equations (^) and 
( |101|) whereas if there is no wall on the axis of symmetry AE% > and AE\ < 0. However 
for low values of no where no and n are close in value to no these conclusions may not hold 
and each case must be treated independently. In particular for (23) + (2) =^ (223) which 
corresponds to the case AE 3 with no wall on the axis of symmetry and for (2 2 3) + (2) =^ (2 3 3) 
which corresponds to the case AE\ with a wall on the axis of symmetry the energy differences 
are negative. 

Hence 0(1/ D 5 ) the phase sequence is 
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(oo) : (4) : (34) : (3); (2333); (233); (23); (23223); (223); (2223) : (2) (127) 

where : denotes a stable boundary and ; a boundary which may be split at higher orders of 
the expansion. 
p = 10 

For p = 10, ao = which means that the leading order terms in the energy differences are 
zero and the series analysis breaks down. Numerically we have been able to show that the 
only stable phases are (2 fe 3) appearing between (2) and (3). The existence of these phases 
has been checked for k < 5 by expanding the energies 0(1/D 6 ). 
p=ll 

A calculation 0(1/ D) shows that the (2) : (oo) boundary is stable; no new phases appear 
near xq and the transition is first order. 



IV. DISCUSSION 

Models with competing interactions may have ground states which include special points, 
so-called multiphase pointsi, where the ground state is infinitely degenerate. At these points 
one might expect that small perturbations can have a drastic effect. One possible such per- 
turbation is to allow the spins to soften, others are temperature^ or quantum fluctuations^!. 

The aim of this paper has been to study the first of these possibilities. We have considered 
continuous spin models with a p-fold spin anisotropy D which in the D — > oo limit exhibit 
a multiphase point. An expansion in 1/D is described which allows us to calculate the form 
of the phase diagram near the multiphase points as the spins soften. 

The first model we considered was the chiral XY model with p-fold spin anisotropy, 
which for large D can be thought of as a soft chiral clock model. This model provided a 
useful illustation of the technique. We showed in agreement with Chou and Griffiths^ that 
for p > 3 all phases formed by combining adjacent structures are stable near the multiphase 
point and obtained leading order expressions for the widths of the phases. 
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We then described results for the more complicated situation of the XY model with first- 
and second-neighbour interactions for large D. Here the situation is very complex with 
the behaviour near the multiphase point being strongly dependent on the value of p. For 
6 < p < 10 infinite sequences of phases are stable, but their quantitative form is different 
for different p. For p > 11 the phase boundary emanating from the multiphase point is first 
order. 

We note that not only does the symmetry of the anisotropy have a non-universal effect 
on the nature of the phase diagram but also that the physical form of the perturbation is 
important. For example, for the ANNNI model itself temperature leads to a phase sequence 
[2 fc 3], k = 0, 1, 2 . . quantum fluctuations result in a phase sequence [k], k = 2,3,4.. E 
whereas when the spins soften there is a single first-order transition. 

The models descibed here have a complex mathematical structure but also interesting 
applications to real systems. For example rare-earth magnetism has been modelled by an 
XY model with competing interactions and 6-fold spin anisotropy^. Modulated structures 
in UNi 2 Si2 have been described using a model expected to behave in a similar way to those 
considered here0. 

Finally we comment on some further possible uses for the technique introduced in this 
paper. Bassler, Sasaki and have descibed an upsilon point, a checkerboard struc- 

ture of long-period phases which is in some sense a two-dimensional version of the sequences 
of phases we have been concerned with here. SasakiH found some evidence for the existence 
of such a point in a spin model by performing an expansion in 1/D to order 1/D 2 . An 
expansion taken to all orders will give firmer proof of where such points can be found0. A 
similar technique can be used to investigate interface unbinding where spin softening can 
allow an interface to unbind from a surface through a series of layering transitions^. 
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APPENDIX A 



We aim here to show how n defined by equations ( j31"D and ( pEUD is related to n[ T ] the 
length of the final phase in the process [a] + \j3\ =>- [7]. 

Consider two odd phases [a] and [j3]. Let no = no ([a/3]) for [a] + \j3\ =>- [a (3\. Using the 
labelling scheme defined in Section |11 G (j) 

[a] = ({«}, -{a}, 0) (128) 

where we use {a} as shorthand for a 2 . . . a( m _i)/ 2 } (see equationfJT^)). Similarly 

[/?] = ({/?},-{/?},0) (129) 

The resulting even phase is 

[a(3] = ({a}, -{a}, 0, {/?}, -{/?}, 0). (130) 

Consider now the process [af3] + \J3] =^ [a 2 /3]. Now an even and odd state are combined 
and therefore the labelling scheme defined in Section [11 C| (ii) is appropriate. From ( |129[ ) and 



[a(3] = (-{{1},0,{a},-{a},0,{{3}), 

\0\ = (-{P}A{P})- (131) 



By inspection of (|131|) it is apparent that 

n ([a(3 2 }) = n ([a{3}) + (n m + l)/2, (132) 

that is adding an odd state [0\ to an even state increases n by (n\p\ + l)/2. 
A similar argument shows that for [a] odd and [/3] even 

n ([a/3 2 ])=n ([a/3]) + (n [/3] )/2, (133) 
n ([a 2 f3}) = n ([a(3}) + (n [a] - l)/2, (134) 
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that is adding an even (odd) state [f3] ([a]) to an odd state [otf3\ increases n by 
(("[«] - l)/2). 

It is not hard to check that the conditions ( |132| ), (|133|) , and ( 134 ) are consistent with 



odd 



n = {n[y\ — (2m + 2)}/2 ri[ 7 ] even 



(135) 



for any integer m. By inspection n ([23]) = 2. Therefore m = and 

"n[ 7 ] + r 



n 



- 1 



(136) 



where [n] is the integer part of n. 
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FIGURE CAPTIONS 



Figure 1: Ground states of the p-state clock model with ferromagnetic first-neighbour inter- 
actions J\ and antiferromagnetic second- neighbour interactions J 2- A vertical line represents 
a first-order boundary and a star a multiphase point. 

Figure 2: Phase diagram of the soft clock model with competing interactions near the 
multiphase point x^ 1 showing the phases that are stabilised by terms in the energy 0(1/ D) 
(labelled above the line) and 0{1/D 2 ) (labelled below the line). A vertical line represents 
a first-order boundary and a star a multiphase line at which an infinite number of phases 
remain degenerate to this order. 
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